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Abstract 

We analyze dynamical instability of non-static reflection axial stel¬ 
lar structure by taking into account generalized Euler’s equation in 
metric f{R) gravity. Such an equation is obtained by contracting 
Bianchi identities of usual anisotropic and effective stress-energy ten¬ 
sors, which after using radial perturbation technique gives modified 
collapse equation. In the realm of i? -|- ei?” gravity model, we investi¬ 
gate instability constraints at Newtonian and post-Newtonian approx¬ 
imations. We find that instability of meridional axial self-gravitating 
system depends upon static profile of structure coefficients while f{R) 
extra curvature terms induce stability to the evolving celestial body. 
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1 Introduction 

General relativity (GR) is believed as a remarkable effort in mathematical 
physics to analyze gravitational effects of stellar relativistic interiors. Sev¬ 
eral interesting consequences coming from cosmic microwave background, 
observational ingredients of Supernovae la and its cross-juxtaposition with 
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foreground stellar galactic distributions n-ia have made revolution thereby 
introducing new research window. In this realm, many astrophysicists found 
GR modihcations as helpful to explore unknown aspects of cosmic gravita¬ 
tional dynamics. The f{R) gravity [1] is among extended gravity theories 
obtained by replacing Ricci invariant with its generic function f{R) in the 
Einstein-Hilbert action. 

Anisotropic effects are leading paradigms in addressing the evolution¬ 
ary mechanisms of celestial imploding models. Herrera and Santos [5] re¬ 
viewed contributions of locally anisotropic fluid arrangements on the dynam¬ 
ical phases of collapsing shear and shear-free compact objects. Di Frisco et al. 
[B] investigated dynamical properties of anisotropic spherical matter distribu¬ 
tion and found that little fluctuations of pressure anisotropy lead to system 
cracking. Sharif and his collaborators [7] analyzed effects of anisotropy on 
the dynamical properties of spherical as well as non-spherical dense rela¬ 
tivistic distributions and found much complicated system phases due to the 
presence of anisotropy. Sunzu et al. [8] studied analytical models of spheri¬ 
cal anisotropic interiors and found that anisotropic effects provide a broader 
platform to discuss various forms of stellar relativistic systems. Recently, we 
[9] explored the dynamical features of anisotropic relativistic interiors. 

The spinning stellar distributions indicate direct relevance of anisotropy 
with gravitational evolution in which gravitational radiations cause vortic- 
ity within observer congruences. Vorticity represents rotation of neighboring 
fluid about an observer moving with relativistic matter distributions relative 
to an inertial frame. Herrera et al. [10] argued that such vorticity seeds 
from the existence of super-energy flow which may have direct relevance 
with super-Poynting vector. Bonnor m found electromagnetic energy flow 
in a relativistic compact distribution by formulating a relationship between 
super-Poynting vector and vorticity. Korunur et al. [T^ calculated various 
kinematical variables like angular momentum, energy and momentum of mat¬ 
ter conhgurations associated with an axially symmetric scalar field. Li [13] 
explored superradiant instability of rotating compact relativistic objects in 
higher dimensional theory and found unstable conhgurations against scalar 
held perturbations. Recently, Herrera et al. [TT] presented a formal analysis 
of gravitational radiations within anisotropic non-static rehection axial sym¬ 
metric source and existence of super-energy how linked with matter vorticity. 

Stability analysis of self-gravitating stellar systems in GR as well as mod- 
ihed gravity have attracted many researchers for last few years. The study of 
diherent collapsing celestial models with extra degrees of freedom has great 
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significance to explore late-time cosmological evolution. Chandrasekhar [T5] 
discussed instability constraints for spherical symmetric relativistic geometry 
coupled with ideal matter conhgurations using ratio of specihc heats known 
as stiffness parameter, hi. Herrera et al. |T6] investigated stability regions 
for radiating collapsing stellar objects and concluded that dissipation vec¬ 
tor tends to move the systems towards stable conhgurations. Chan et al. 
El studied remarkable effects of shearing viscosity and anisotropy on the 
instability constraints at Newtonian (N) and post-Newtonian (pN) eras. 

Cai [18] discussed dynamical properties and structure formation of dense 
matter relativistic conhgurations in modihed gravity by assigning zero, nega¬ 
tive, and positive values of constant curvature. Bamba et al. [12] performed 
dynamical analysis of a collapsing relativistic stellar system and claimed that 
invoking of i?"(l < a < 2) corrections could helps to present a viable and sin¬ 
gularity free model. Myung et al. [20] performed stability analysis of spheri- 
calstellar structure with constant Ricci invariant background in metric f{R) 
gravity via perturbation scheme and noticed relatively stable distributions 
under specihc constraints. Moon et al. [2T] extended these consequences 
with negative cosmological constant environment and calculated limits for 
the stability of relativistic systems. 

Capozziello et al. [22] explored dynamical evolution of relativistic collaps¬ 
ing spherical interior in f{R) gravity by evaluating extended form of Poisson 
and Boltzmann equations. De Laurentis and Capozziello [23] discussed in¬ 
stability issue of stellar interior at N approximation with f{R) extra degrees 
of freedom and also studied axisymmetric black hole models. Astashenok et 
al. [23] investigated evolution of self-gravitating systems and found relatively 
more massive and supergiant dense conhgurations due to f{R) gravity cor¬ 
rections. Farinelli et al. [23] discussed dynamical properties of stellar systems 
in the presence of f{R) corrections and found that higher degree terms tend 
to mollify collapsing process. Sharif and his collaborators [26] 127] studied 
instability constraints for restricted class of axially symmetric spacetime by 
means of adiabatic index/stihness parameter. 

The present paper aims to extend our previous work m of stability 
analysis by taking rehection ehects in non-static axial symmetric anisotropic 
source with ei?" extra degrees of freedom. In the present paper, we de¬ 
velop instability regions for anisotropic meridional axisymmetric source with 
R -|- eR^ background. The inclusion of eR^ correction in our analysis seeds 
from the fact that this correspond to the various eras of the cosmic history 
thereby helping to explain the gravitational dynamics during inhationary as 


3 


well as late-time accelerating regimes. Furthermore, the addition of merid¬ 
ional effects in stellar system causes a flow of gravitational energy due to 
existence of vorticity tensor in the analysis. 

The paper has the following format. Section 2 deals with kinematical 
formulations of comoving meridional axial symmetric geometry coupled with 
anisotropic matter conhgurations. The meridional effects in stellar system 
causes a flow of gravitational energy due to existence of vorticity tensor. 
We present f{R) dark source components and set of dynamical equations 
with reflection axial degrees of freedom. In section 3, we discuss viable f{R) 
model and use perturbation method to develop collapse equation. Section 4 
explores instability constraints. Finally, we summarize our results in the last 
section. 


2 Anisotropic Source and Field Equations 


The extended conhguration of Einstein-Hilbert action is 

Sf{R) = + Sm, ( 1 ) 

where k, f{R), Sm, Ta0 are coupling constant, matter action, a non-linear 
Ricci curvature function and usual stress energy tensor, respectively. The 
variation of above action with respect to gap provides the held equations 

IrRolP — pIr — dap (^2^ ~ = I'^Tap, ( 2 ) 

where □, Va, fa are D’Alembert, covariant derivative and ^ operators, 
respectively. Equation ([2]) can be written in terms of Einstein tensor as 


K (^) 

Gap = -^{Tap Tap), (3) 

JR 


where 

S = ^ {f (^ - fR^ 9ap - Ofngap + VaV^/^l , (4) 

is the stress energy tensor indicating f{R) contribution in the dynamics of 
relativistic systems. We take axially symmetric metric characterizing rehec- 
tion effects [H] 

r, 6)dR + 2L{t, r, 6)dtd6 + B'^{t, r, 6){dr‘^+r'^d6'^)+C'^{t, r, 6)d(j)^, 

( 5 ) 
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with locally anisotropic fluid conflguration 

Tap = (h + P)Vayp + PgaP + TPl){KaKp ~ 211//) 

X {NaNp — -hap) + T\-KN{KaNp + KpNa)-, (6) 

where p, P, 11/, 11//, 11/^/, and hap are the fluid energy density, pressure, 
anisotropic scalars and projection tensor, respectively. The matter four ve¬ 
locity, Vai and spacelike vectors Sa, Ka and Na in comoving coordinates 
are 


v; = -.4^" + !^sl, = cy, a'„ = bsI n,, = ^sl (7) 


where A = [ABr)^ -|- L^, which obey the following constraints 


K^Na = K^Sa = S^Na = VaK^ = V^Na = V^Sa = 0 , 

KaK^ = NaN° = SaS^ = -V'^Va = 1. 


The fluid pressure and its anisotropic scalars can be expressed alternatively 
in terms of projection tensor and spacelike vectors, respectively as 

P = U^^Tap, n/ = (2K^K^ - - N^N^)Tap, Ukn = K^N^Tap, 

Ujj = - K'^K^ - S°S^)Tap. 

The non-zero components of effective stress energy tensor (jl]) are 

1/i + lTi Ai + Fi X 3 + F 3 0 

Ai + Fi V 2 + W 2 X 2 + Y 2 0 

A 3 + 1^3 A 2 + F 2 V^3 + W^3 0 ’ 

0 0 0 14 + 

where dark source /(P) terms Vj, IT* and Xj, Yj are diagonal and non¬ 
diagonal components of effective energy-momentum tensor (jlj), respectively, 
in which Wi and Yj incorporate axial reflection effects with /(P) extra degrees 
of freedom. By choosing Xj and Yj equal to zero along with A —)■ 
higher curvature terms of restricted axisymmetric metric can be found. How¬ 
ever, the inclusion of these terms along with anisotropic in usual stress energy 



5 





tensor ensure the propagation of gravitational radiations in the environment 

izg. 

The kinematical quantity controlling local spinning motion of anisotropic 
matter distributions is the vorticity tensor which for meridional axially sym¬ 
metric metric can be expressed in terms of Ka and Na as 




where 




L fV 2kl'\ 
2By/K A J 


(9) 


is known as vorticity scalar. Here prime stand for There exists only one 
independent non-zero vorticity component along rO direction. The existence 
of vorticity scalar is directly related to the existence of reflection effects of 
axisymmetric spacetime as it is controlled by non-diagonal structure coef¬ 
ficient, L. Thus if one takes H = 0 over the dynamical evolution of axial 
anisotropic spacetime, this imparts null value to non-diagonal scale factor 
whose dynamics has already been discussed |Z7|. 

In order to evaluate dynamical evolution equations for axially symmetric 
relativistic celestial body with f{R) background, we consider 


(D) 

|T"^ + r“^|.^ = 0, 
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( 12 ) 


where Dq, Di and D 2 are f{R) corrections given in Appendix A. Here over 
dot and subscript 6 stand for ^ and respectively. The second of the 
above equations is known as generalized Euler equation. 


3 f{R) Model and Perturbation Scheme 

Many inflationary models in the early universe are established on scalar fields 
coming in supergravity and superstring theories. The hrst model of inflation 
was suggested by Starobinsky which deals with conformal anomaly in quan¬ 
tum gravity [29] given by [30] 

f{R) = R + eR!', (13) 
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where n can be positive or negative integer. This model explains the present 
universe acceleration due to the presence of dark energy and can serve as 
power-law inflation, i.e. exponential expansion and ordinary inflation incor¬ 
porating minimally coupled scalar held. Here e ~ jyj 2 n -2 > 0 for n > 0 and 
M has the mass dimensions. Since f{R) gravity can be used as an alter¬ 
native for dark matter |3I] in addition to dark energy at cluster as well as 
stellar scales, so this model with n = 2 was claimed both as dark matter 
model with e = [32] and as dark energy. The value of M is chosen to be 

2.7 X along with e < 2.3 x for dark matter cosmology 

[33] . All GR solutions can be found by taking limit f{R) —)■ R. 

Here, we use perturbation method [SHU] to explore modihed collapse 
equation for meridional axially symmetric anisotropic geometry. For very 
small values of perturbation parameter a with 0 < a <C 1, we take effects 
up to 0{a). We hrst suppose that the system is in hydrostatic equilibrium 
at f = 0, however on departing from this state, the system depends upon 
the same time dependence factor T{t) on all its structure coefficients. The 
structure and matter variables can be perturbed as follows 

S{t,r,e) = So{r,e) + aT{t)s{r,e), (14) 

M{t,r,9) = Mo{r,9) + am{t,r,9), (15) 

where S represents perturbation method applicable on structural co-efficients 
of Eq.(|5]), i.e.. A, B, C, L and on Ricci scalar, R which after perturbation 
denotes s as a, b, c, I and e, respectively. Equation fITBD indicates pertur¬ 
bation method of matter variables (these matter variables are taken from 
Eq.((6])). Thus the allocation of Ai will be /r, P, Hq, a = 1,2,3 and the 
corresponding perturbed quantities will be represented by placing bar over 
that. However, the perturbation technique for f{R) model is given as follows 

r) = [Ro{r) + eRo)] -F aT{t)e{r) [l - enR^-^] , (16) 

fnityr) = 1 + enRQ~^ + aT{t)e{r)ne{n — 1 )Rq~‘^ , (17) 

where Rq represents static distribution of Ricci scalar. Using Eqs. ffTT|) - ffT7|) . 
the first of dynamical equations satisfies trivially, while rest of dynamical 
equations flTTD and (IT^ at f = 0 give 

-Pq + g (2njo + n^jo) + 
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The static f{R) contribution of second and third conservation equations 
are denoted by -D 25 and -D 35 , respectively and can be calculated very easily 
from Eqs. flA2l) and flA3ll after using perturbation method. Using Eqs. (ll4p - 
m, the non-static perturbed axial dynamical equation flTOl) will take the 
form 
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where D 3 represents f{R) corrections which can be obtained from expressions 
g{t,r,9) and h{t,r,9) given in Appendix A. Substituting Eq. flT^ in Eq.(| 8 ]) 
and then employing perturbation method, one can obtain f{R) dynamical 
quantities, V), Wi, Xj, Yj whose values upon substitution in Eqs. flAdp and 
flASp yield -D 3 such that g{t, r, 9) + h(t, r, 9) = D^T. Integration of the above 
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equation gives 
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Now, we evaluate t9 component of metric f{R) field equations ([2]) and 
then using perturbation scheme along with some manipulations, it follows 
that 


PiT + + QsT — 0, 


( 21 ) 


where quantities pj contain combinations of meridional axial geometric func¬ 
tions as well as i? -|- ei?" corrections, depending upon r and 9 coordinates 
and are assumed positive. More specifically, these quantities incorporate 
non-perturbed as well as perturbed terms. There exist oscillating as well as 
non-oscillating solutions of the above equation which represent unstable as 
well as stable models of evolving relativistic stellar systems, respectively. We 
conhne ourselves to obtain solutions for collapsing relativistic system. Thus 
we limit our perturbation parameters, a, 6, c, e and I to be positive defi¬ 
nite quantities for which we obtain > 0. In this context, the solution of 
Eq-dZID is given by 
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where T, and $ are mentioned in Appendix A. The qnantity controlling the 
reflection degrees of freedom along with f{R) corrections of an axisymmetric 
celestial body is T. However, the expression $ incorporates gravitational 
contribntion dne to f{R) gravity, while C, is the remaining part of non-static 
pertnrbed generalized Enler eqnation holding usnal Einstein gravity effects. 

In view of second law of thermodynamics, we can link pertnrbed anisotropic 
qnantities with energy density by an eqnation of state as [5Tj 
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where Ti is a flnid stiffness parameter also known as adiabatic index. This 
measnres pressnre variations of matter conflgnrations with respect to energy 
density. In onr analysis, Ti will be treated as a constant identity. Using 
Eqs. (l2n|l and (Elj), we have 
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( 25 ) 


= (j£m + a(™”n„) + sIm+SL) ' equation is known as 

collapse equation of axisymmetric stellar objects characterizing meridional 
and /(-R) extra order degrees of freedom. 

4 Instability Regions 

Now we proceed to calculate constraints at which meridional axial symmet¬ 
ric stellar systems undergo instability window at both N and pN eras with 
f(R) background. We also examine the role of stiffness parameter hi in 
this scenario. We also reduce our results to previously known limiting cases. 
The formulation of instability constraints should be compatible with Tolman- 
Oppenheimer-Volkoff (TOV) equation. Such type of equation constrains the 
relativistic stellar structure coupled with matter distribution at the phase 
of static gravitational equilibrium. In this respect, Barausse et al. [35] in¬ 
vestigated hydrostatic equilibrium phases of relativistic models by obtaining 
modihed version of TOV equation in f{R) gravity. Recently, Astashenok 
et al. [3n| calculated extended version of TOV equation with equation of 
state in the realm of cubic as well as quadratic corrections and found that 
such an equation can be used to describe viable models of compact objects. 
Here, we formulate TOV equation that will help us to obtain some limits on 
fluid-energy density and its derivatives to avoid curvature divergence at the 
stellar boundary. The 11 and 22 components of metric f{R) held equations, 
respectively, provide 
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The corresponding Misner-Sharp mass function [S^ takes the form 
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Using Eg. (1301) . second and third laws of conservation of usual energy-momentum 
tensor as well as Eqs. (12611 and ([27]), we obtain TOV equations 
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where ijjm is calculated by taking an assumption that reflection effects are far 
lesser than that produced by other scale factors in the evolution of axisym- 
metric system. In order to examine the contributions of f'^ and f^g across 
the meridional non-static axial relativistic object, we multiply both sides of 
the above equations with After some manipulations, this yields couple 
of quadratic equations in and f^g whose solutions become 
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where Ai and A 2 are discriminants of f'p^ and fRg quadratic equations, C = 


9^ir^A^■^|J^-9r^^A^'^P^^A{9P+2{nII+2UI)} 

9A2 


Cl = 2{Un+2Ujy^, while D = ±1. 
We shall take Aq = 1 — ip, Bq = 1 + (p with ip = ^ for pN epochs, therefore 


Over the surface of axial reflection relativistic star object, Eqs. fIST]) and (1321) 
yield 


fy ^ W,{D-1) ^ W2{D-1) 

fR iscy^mA ' fR 2 ’ ^ ^ 

where Wi = 144(i/>mO rCxp'^) and W 2 = + 

2LC'il)me - 2LCgiJm) - A^B‘^r^K{BCg - CBg) + r^CA‘^{r^A^B^ - L^). It can 
analyzed from Eq. fl33l) that on setting D = —1, one can get specific forms 
of 7 and 7 i from Eqs. fl2^ and (l2^ which will make ^ and ^ approach to 
00 with (r, 6 ) —)■ (r“, 6 ~), while finite value of ^ and ^ can be achieved for 
(r, 0) —)■ For physically viable stellar model, we take D = 1 which 

yields fRg = 0 = /jj for (r, 0) {r~,9~). This reinforces the continuity of 

fne, I'r as well as A' over the surface of axial stellar structure with reflection 
degrees of freedom. 
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4.1 Newtonian Approximation 


In order to evaluate instability conditions at N regime, we take Aq = Bq = 1 
and assume anisotropic pressure to be less than zero which is the criterion 
for collapsing celestial body. We also take conhgurations of initial perturbed 
structural coefficients to be Cq = Lq = r. Consequently, the collapse equation 
turns out to be 


9 


n 


Km,e 


/2c 


+ jr-A.vr. - r. ^ ^ 


— + 36 + 


3 

8r 


^ + 36 + ^ ) - 26 [P' 


-m'j. + U'jj,) 


+ 


2rV2 V 


( 2b -^ + -Xsoat + —Xioat + P + + Cj 


2r 


where subscript J\f indicates the evaluation of term at N regime. We assume 
that all terms on both sides of the above equation are positive. The instability 
constraint for meridional axisymmetric fluid configurations is given by 


ri< 


(y + 36 + ^) - 26 [Pq + 1(211^0 + Hj/q)] + 0i + Cv 

(f + 3^ + 


(34) 


where 0i = fXgoAr + ^Xioat + ^ (26 - ^ + Tat + <hAt + Cjv- The sys¬ 
tem would be in complete hydrostatic equilibrium, if (during evolution) it 
can take a value equal to the right hand side of the above expression. How¬ 
ever, on satisfying the above inequality, the system will move in the unstable 
phase. This constraint is being mentioned through Ti parameter thereby 
emphasizing the importance of matter stiffness factor in our investigation. 


4.2 Post-Newtonian Approximation 


Here, we take axial structural coefficients for pN eras and consider our out¬ 
comes upto 0{ip). Using these relations in Eq. fl25|) . one can have modified 
collapse equation at pN limit. This leads to instability inequality through 
stiffness parameter 


Ti < 


r2(l - 40){(p' + i(l - (p)(l - ^)e}XpN + (1 - 2(p)^((p' + |:) -|- Ci 


(1 - 2(p)0W + (1 - 2(p)0pAr[^ + 2(r ~ + C 2 


,(35) 


where 


Cl — Axme — 26(1 — 3ip) 


Tq + g(2n/o + n'jjo) 


+ 


(l + 2 (p) 

rx /2 


(36 — a + 2aip 
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31 


—3b(p — — j + —XiopAT + 
r^{l — (p) HknoXpN 

^2 — ■ 


6(1 - V?) 


X. 


30pN 


+ '^pN + ^pN + CpNi 


[6(1 + (^)e(l ~ 9^) + (1 + ^)o{^ + ^) + 2(1 + 


2\/2 /io + Hi^Tvo 

X{(1 - »,)„(1 +»,) + (1 + ^)„(1 - m (iki^) , 

zyz V/^o + e 


the subscript pN represents effects of quantities at pN era. The quantity 
Tp 7 v describes the reflection effects of non-static axial celestial body about 
its symmetry axis at pN approximations. It is worth mentioning here that 
these constraints coincide with 123 in the limit L —)■ 0 for n = 2. 


5 Instability of Realistic Star Object 

Perturbations of stars and black holes have been one of the main topics of 
relativistic astrophysics for the last few decades. The description of such 
stellar objects has recently attracted various researchers [38]. The stability 
analysis of general relativistic star process is an important but challenging 
endeavor. In such study, the spherical symmetric matter conhguration is 
an exemplary one. Numerous realistic objects like globular clusters, galactic 
bulges and dark matter haloes can be considered as being roughly spherical 
geometry. For better understanding of cosmic censorship hypothesis and 
hoop conjecture, it is necessary to throw light on non-spherical collapse. The 
physical interest in studying non-spherical symmetries is associated with the 
fact that post-shocked clouds are left at the verge of gravitational collapse 
forming cylinders or plates at scales of galaxy formation and at scales of 
stellar formation in galaxy. For instance, cylindrical distributions are closely 
related with the problem of fragmentation of prestellar clouds [39] . 

We take into account a specific conhgurations of non-static axial space- 
time. The main purpose is to study instability regimes of axially symmetric 
realistic objects that are involved in the emission of gravitational radiation 
due to meridional degrees of freedom. For this purpose, we assume coupling 
of system with anisotropic fluid distribution whose energy-momentum tensor 
is mentioned in Eq.(|5]). Having arrived at this point, the relevant question is: 
does an ideal (or non anisotropic) matter configuration produce gravitational 
radiations? 

To answer such a burning issue, let us recall that in the seminal paper of 
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Bondi about the emission of gravitational radiation (section 6 of [ID]), it is 
mentioned that for relativistic dust cloud as well as dissipation-less case of 
an ideal matter distribution, the relativistic system cannot be anticipated to 
radiate (gravitationally). This is due to the reversible feature of equation of 
state as emission of radiation is an irreversible phenomenon. This happens 
once when absorption is considered (and/or Sommerfeld type constraints), 
which prevents inflow of waves. This implies that an entropy generator pa¬ 
rameter/factor must be present in the discussion of relativistic source. How¬ 
ever, such type of factor is not present in an ideal fluid and in a collisionless 
dust cloud. In particular, the irreversibility of gravitational wave emissions 
must be taken in equation of state with the help of an entropy increasing 
(dissipative) parameter. In this scenario, Herrera et al. [TT] described a close 
relationship between vorticity and gravitational radiations. 

We consider the evolution of non-static axisymmetric self-gravitating sys¬ 
tem in f{R) gravity and assume that it is in hydrostatic equilibrium at large 
past time. Now we want to analyze that when the phase of equilibrium is 
disturbed, what happens? Will this perturbation be relaxed (stable state) or 
will it grow (unstable state). In this respect, one needs to take into account 
couple of following instabilities 

1 . dynamical stability: what happens, if stellar hydrostatic phase is per¬ 
turbed? 

2 . secular (thermal) stability: what happens, when the state of thermal 
equilibrium is perturbed? 

Since our system is coupled with anisotropic matter conhgurations without 
heat flux, therefore we shall not discuss the second case and conhne our¬ 
selves over dynamical instability of relativistic origin. It is seen that under 
hydrostatic phase, the stability criterion is achieved by making linearized 
held as well as conservation equations against radial perturbation (14)-(17). 
It is remarked that during evolution, the realistic object moves via several 
evolutionary patterns determined by instability/stability degrees of freedom. 
This suggests that the relativistic systems can be stable at one instance but 
not at the other. Thus one needs to cope with the dynamical evolution 
of self-gravitating systems by calculating instability regions at N as well as 
pN regimes. Such epochs have vital role in the discussion of gravitational 
collapse of compact objects. 
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The phenomenon of celestial collapse occurs when the state of hydrostatic 
equilibrium of a stellar object is disturbed. In celestial body, nuclear hssion 
reactions occur that start from hydrogen atoms and produce further com¬ 
plex elements until nuclear reactions chain stops with iron. These reactions 
increase the pressure exerted by gas particles which counterbalance the grav¬ 
itational attraction and prevents the star from collapsing. However, with the 
passage of time, nuclear reactions decrease as fuel burns out. Consequently, 
the necessary pressure becomes insufficient for a collapsing body to be stable. 
At this point, the gravitational force begins to pull matter towards the center 
of a body and thus collapse initiates. A celestial body that has exhausted 
all its nuclear fuel, can give birth to three possible compact objects (white 
dwarfs, neutron stars and black holes) on the basis of the initial mass of the 
collapsing body. 

It is well-known that, in the scenario of Newtonian regime, the instabil¬ 
ity of spherical self-gravitating systems depends purely on the mean value 
of stiffness parameter, Ti |12] which is the ratio of fractional Lagrangian 
variations between pressure and energy density experienced by matter con- 
hgurations following the motion. However, in GR, the stability relies not 
only on the average value of Ti but also on the star radius. However, in the 
study of non-static axial reflection system in modihed gravity, the situation is 
quite different. (It is worthy to stress that we have assumed Ti as a constant 
entity throughout the analysis). The most important consequence of our 
study is that, apart from affirming GR results, Ti controls emission of grav¬ 
itational radiations along with f{R) extra degrees of freedom. The emission 
of gravitational radiations causes the loss of both energy as well as angular 
momentum which increases the instability of the meridional axisymmetric 
object. 

More specihcally, following the results of Chandrasekhar [15], we deduce 
that if the anisotropic matter distribution attains stiffness equal to the right 
hand side of expressions fIMD and fl35H . the system enters into the window of 
hydrostatic equilibrium at N and pN regimes. Further, if the stiffness of fluid 
increases in such a way that the fractional value given at the right hand side 
of fl5T|) and fl35p becomes a smaller one, then system enters into the stable 
conhgurations at both N and pN approximations, respectively. Dosopoulou 
et al. |l3] explored the contribution of magnetic helds in the emergence and 
existence of vorticity. This strongly suggests that invoking of magnetic helds 
in the study of stability of gravitationally radiating sources deserves attention 
for future work. 
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6 Conclusions 


It is well-known that the most general non-static axial geometry incorporates 
reflection (meridional) and rotation effects coming ont from non-diagonal 
dtd6 and dtdcj) metric coefficients. In order to dealt analytically with in¬ 
stability constraints of axially symmetric spacetime, several attempts have 
been made by taking restricted class of axial geometry. In this paper, we 
have studied stability analysis of meridional axial stellar structure with f{R) 
background. We are observing investigation in a metric f{R) gravity which 
give rise to non-linear fourth order held equations. We have formulated the 
collapse equation by using perturbation scheme in the generalized Euler equa¬ 
tion. We assume complete hydrostatic equilibrium of axial stellar structure 
at large past time, i.e., T(—oo) = 0. 

We have developed instability constraints at N and pN epochs through 
stiffness parameter, Ei using collapse equation. It is found that axial stellar 
structure would be unstable until it obeys relation fl34p at N regime while 
relation 0351) at pN era. Breaching of these inequalities will eventually move 
the system towards stable window. These constraints depend upon adiabatic 
index, static combinations of anisotropy, energy density and dark source cor¬ 
rections due to R + eR^ model. It is seen that dark source corrections tend to 
stabilize structure formation phenomenon due to its non-attractive behavior 
while the presence of non-diagonal terms in instability ranges indicate oc¬ 
currence of gravitational radiations which correspond to flow of super-energy 

[H]. 

We have found non-vanishing component of vorticity tensor which cor¬ 
responds to non-static meridional axial structure coefficient. The inclusion 
of non-diagonal scale factor in the stability analysis leads to interesting phe¬ 
nomenon of gravitational radiations for eR!^ corrections. These extra-order 
f{R) corrections affect the passive gravitational mass which in turn affect the 
rate of stellar collapse. We have developed instability constraints fl3T)) and 
fl35|l with weak held and pN approximations, respectively. These constraints 
can be applied to axisymmetric self-gravitating system with reflection degrees 
around symmetry axis at some particular cosmic epochs depending upon the 
chosen values of n. We can categorize different eras of cosmic dynamics 
associated with ei?” as follows. 

• For n = 2, the instability constraints for specific model of the type 
R -|- eR^ can be obtained. The existence of R^ correction in the held 
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equations can be helpful to explain inflationary mechanism of cosmos. 
The term aR^ represents accelerated expansion of the universe. This 
model is compatible with temperature anisotropies noticed in cosmic 
microwave background radiations [3l] and hence viable for inflationary 
scalar held models. 

• The choice n = 3 favors to host signihcant massive compact objects 
coming out from cubic f{R) higher curvature terms |2Tj. This pro¬ 
vides realistic signature of the presence of more massive and huge self- 
gravitating stellar systems which have direct correspondence with the 
observational cosmology. 

• This gravitational dynamics at late-time universe era can be obtained 
by substituting u = — 1 in instability constraints at both N and pN 
regimes. It is noticed that gravitational contribution due to negative 
curvature power serves as dark energy thereby supporting current ac¬ 
celerating cosmic epochs [H]. 

• For e = 0, instability constraints for Einstein gravity can be obtained 
at both N and pN eras which describes relatively less stable axial stellar 
structure. 

Finally, we remark that supermassive stellar systems survive more abun¬ 
dantly in extended gravity than in GR as such theories (for instance f{R) 
gravity) are more likely to host huge stars with smaller radii. This leads to 
the existence of more dense relativistic systems which have direct relevance 
with observational gravitational physics. It is interesting to mention here 
that all our results reduce to restricted class of instability analysis m by 
neglecting non-diagonal terms and assuming u = 2 in f{R) model. 
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Appendix A 

The extra f{R) curvature terms for Eqs. fllOI) - (fT^ are 
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